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INTRODUCTION 



Based on many observations in hadron physics spontaneous breaking of chiral symmetry 
and the dynamical generation of quark masses are expected to occur in Quantum Chromo 
Dynamics (QCD). The precise origin of this non-perturbative phenomenon as well as its 
relation to quark confinement are still little understood. Further studies of these issues 
have to build on reliable non-perturbative methods. And, as confinement is expected to be 
correlated with infrared singularities, continuum-based methods will be required in addition 
to Monte Carlo lattice calculations. To this end we note that the Schwinger-Dyson equations 
(DSEs) of QCD can address directly the infrared region of momentum. 

The DSEs for the propagators of QCD form a coupled system of equations. In Landau 
gauge these have been investigated in two lines of research. On the one hand, the DSEs of 
pure Yang-Mills theory have been explored with the aim to reveal the infrared behavior of the 
ghost and gluon propagator and their relation to gluon confinement (see e.g. the review |]J 
and references therein). On the other hand, the quark DSEs have been studied extensively 
on the purpose of model building. These models have been used in the frameworks of Bethe- 
Salpeter equations and finite temperature field theory to describe hadronic properties and 
reactions in a semi-phenomenological way (see e.g. the review |2j and references therein). 

In this article numerical results for the coupled set of quark, gluon and ghost propaga- 
tors, including the backreaction of the quarks on the ghosts and gluons, are presented for 
the first time. In the quark DSE we will study the mechanism of dynamical chiral symmetry 
breaking by which physical quark masses are generated even though the bare quark masses 
in the Lagrangian are zero. This is a genuine non-perturbative effect as it is well known that 
for vanishing bare masses the renormalized masses remain zero at each order in perturbation 
theory. In addition to the phenomenon of mass generation we are interested in quark confine- 
ment. Single quark states have non- vanishing color charge and are therefore not contained 
in the physical part of the state space of QCD. This physical subspace supports a positive 
(semi-)definite metric whereas the remaining state space of QCD contains negative norm 
states as well. Consequently, negative norm contributions to the quark propagator would 
provide evidence for quark confinement. 

This paper is organized as follows: In the next section suitable ansdtze for the quark- 
gluon vertex will be constructed such that the DSE for the quark propagator guarantees 
the realization of two important consistency conditions. These are (i) the independence 
of the dynamically generated quark mass function from the renormalization point and (ii) 
the correct asymptotic behavior at large momenta such that the anomalous dimensions of 
dressing and mass functions are correct in one-loop order. Fortunately, the corresponding 
DSEs for the fermions of QED are well studied (a short overview is given e.g. in j^). We will 
dwell on these results and construct non-Abelian generalizations of Abelian vertices, which 
have the desired properties. 

In the following section we present solutions for the quenched system of quark, ghost and 
gluon DSEs, i.e. we neglect the quark-loop in the gluon equation. The Yang-Mills sector is 
hereby treated in a truncation scheme for the ghost and gluon equations that already has 
been employed in ref . |4] . This scheme improves on older ones p, 0| , provides an explicit 
numerical solution for the infrared analysis given in refs. 0, H| , and its results are in almost 
quantitative agreement with corresponding results of lattice calculations for the gluon and 
ghost propagators 0, [HI 11, IHl- The main purpose of calculating the quark propagator 
in quenche d ap proximation is to allow for a comparison with corresponding recent lattice 



results |13l Il4j . We find very good agreement of our results for the quark renormalization 
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function and the momentum- dependent dynamical quark mass with the lattice results if the 
quark-gluon vertex functions constructed in Sect. 2 are employed. We then proceed to the 
unquenched case and incorporate the quark-loop into our truncation scheme for the ghost 
and gluon DSE. We present solutions for the full coupled system of DSEs for the quark, 
ghost and gluon propagators. Compared to the quenched case we will find only moderate 
differences for the number of light flavors Nf < 3. Whereas we are able to demonstrate 
positivity violation in the gluon propagator (and thus gluon confinement) we have not been 
able to draw definite conclusions on the (non-)positivity of the quark propagator. In the last 
section we summarize our results and present some conclusions. 



THE QUARK DYSON-SCHWINGER EQUATION 

The coupled DSEs for the gluon, ghost and quark propagators in Landau gauge QCD 
are shown in in fig. ^ As stated in the introduction we will employ first the quenched 
approximation, i.e. we will neglect the quark-loop in the gluon equation. This will allow us 
to assess the quality of our quark DSE solution by comparing to lattice calculations of the 
quark propagator [l3j performed so far only in quenched approximation. As we furthermore 
employ the gluon propagator obtained in ref. the unknown element to be determined in 
the quark DSE is the dressed quark-gluon vertex. 



Following the conventions and notations of ref. JlJ the renormalized quark DSE with 
appropriate quark wave function and quark-gluon vertex renormalization constants, Z 2 and 
Zip, respectively, reads 



S-\p) = Z 2 S^(p) + g 2 Z lF C F j^l, S(q) T u (q, q - p) D, u (q 



p) . 



The factor Cp = (N 2 — 1)/2N C in front of the integral stems from the color trace of the loop. 
The symbol T u (q, q — p) denotes the full quark-gluon vertex. Suppressing color indices the 
quark and gluon propagators in Landau gauge are given by 

S(v) = 1 = ,1- W^ +M(p2) (2) 

Ki > -ij> A{p 2 ) + B{p 2 ) VP } p 2 + M 2 {p 2 ) ' KJ 

Soip) = \Z~~i ' ( 3 ) 

—ip + mo 

D M =(i„-»f)*f., (4) 

where the quark mass function M is defined as M(p 2 ) = B(p 2 )/A(p 2 ). 

The inverse of the vector self-energy, 1/A(q 2 ), is often denoted as the 'quark wave function 
renormalization'. The bare quark propagator So(p 2 ) contains the unrenormalized quark 
mass m (A 2 ) which depends on the cutoff A of the theory. The bare mass is related to the 
renormalized mass rriR(fi 2 ) via the renormalization constant Z m \ 

m (A 2 ) =Z m (/i 2 ,A 2 )m fi ( /U 2 ). (5) 
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FIG. 1: Diagrammatic representation of the Dyson-Schwinger equations for the gluon, ghost and 
quark propagators. The wiggly, dashed and solid lines represent the propagation of gluons, ghosts 
and quarks, respectively. A filled blob represents a full propagator and a circle indicates a one- 
particle irreducible vertex. 
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Here fi 2 is the (squared) renormalization scale. Note that the propagator functions A, B and 
Z depend on this renormalization scale. Therefore in the following we will use the notation 
A(p 2 ,ii*),B(p2,^) and Z(p 2 ,^). 

The renormalized and the bare vector self-energy, A and A , are related by 

A \p 2 ,A 2 ) = Z 2 ^,A 2 )A-\p 2 ,^). (6) 

In Landau gauge the loop corrections to the vector self energy are finite. Correspondingly 
Z 2 (/i 2 ,A 2 ) stays finite when the cutoff is sent to infinity, and we have < Z 2 (fi 2 ,A 2 ) < 1. 
Furthermore, in Landau gauge the ghost-gluon vertex is not ultraviolet divergent, and we can 
choose Z\ = 1 2|- The Slavnov-Taylor identity for the quark-gluon vertex renormalization 



factor Zip thus simplifies (l6j | 

„ ZiZo Zi ._. 
Z 1F = ^1= * (7) 

Previous studies of the quark DSE in the so-called Abelian approximation (see e.g. the 
review or the recent summary |l7| and references therein) as well as the recent investi- 
gation in ref. assume implicit cancellations between the full quark-gluon vertex and the 
dressed gluon propagator in the integral over the kernel of the quark DSE. Furthermore, in 
the tensor structure of the quark-gluon vertex only a term proportional to 7^ is employed. 

In the current investigation we do not have to rely on implicit cancellations since we 
calculated explicit solutions for the dressed gluon and ghost propagators pj. We will also 
construct explicit non-perturbative ansdtze for the quark-gluon vertex. We note that very 
recently lattice results for the quark-gluon vertex became available 0]. However, at present 
the error bars from such simulations are too large to use the lattice results as guideline in 
the construction of reliable ansdtze for the quark-gluon vertex. To proceed we assume that 
one may approximately factorize the quark-gluon vertex 

r„(<7, k) = Vf el {p, q, k) W- abel (p 2 , q 2 , k 2 ), (8) 

with p and q denoting the quark momenta and k — (q — p) the gluon momentum. The non- 
Abelian factor \y^ abel multiplies an Abelian part V® bel , which carries the tensor structure 
of the vertex. This ansatz is motivated by the aim to respect gauge invariance as much as 
possible on the present level of truncation. 

The Slavnov-Taylor identity (STI) for the quark-gluon vertex is given by 2 



G-\k 2 ) 1 k, r M (g, k) = S-\p) H(q,p) - H(q,p) S'\q), (9) 

with G(k 2 ) being the ghost dressing function and H(q,p) the ghost-quark scattering kernel. 
At present the non-perturbative behavior of the ghost-quark scattering kernel is unknown. 
Therefore we cannot solve this STI explicitly. However, comparing the structure of eq. (JHJ) 
with the corresponding Ward identity of QED, 

tk ll T^ D (q,k) = S- 1 (p)-S-\q), (10) 

we are able to infer some information: Whereas the ghost fields of QED decouple from the 
theory and consequently do not show up in the Ward identity there is an explicit factor of 
G^ik 2 ) on the left hand side of eq. ©• We therefore suspect the quark-gluon vertex of QCD 
to contain an additional factor of G(k 2 ) compared to the fermion-photon vertex of QED. 
Some additional ghost dependent structure seems necessary to account for the ghost-quark 
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scattering kernel on the right hand side of eq. (|9~]l. For simplicity we assume the whole ghost 
dependence of the vertex to be contained in a non-Abelian factor multiplying an Abelian 
tensor structure. For quenched calculations and in the context of angular approximated 
DSEs a similar strategy has already been adopted in refs. 0, 21 j . The Abelian part of 
the vertex, V^ bel , will be adopted in the following from analog ansatze in QED. The Ward 
identity (|1U|) has been solved e.g. in ref. such that kinematical singularities are avoided. 
Furthermore transverse parts of the fermion-photon vertex have been fixed such that mul- 
tiplicative renormalizability in the Abelian fermion DSE is satisfied for all linear covariant 
gauges |23| . The resulting vertex is known as the Curtis- Pennington (CP) vertex. 

The non-Abelian factor \y^ abel is chosen such that the resulting quark propagator fulfills 
two conditions: 

(i) The quark mass function M(p 2 ) should be independent of the renormalization point 

(ii) The anomalous dimension 7 m of the mass function known from perturbation theory 
should be recovered in the ultraviolet. 

In the course of this section we will prove the vertex ansatz 

/ \ ~2d-d/5 

(G(fc 2 ,/i 2 )Z 3 (/i 2 ,A 2 )) 

1 j ^ ' K ' (Z(k 2 ,fi 2 )Z 3 (fi 2 ,A 2 )) d 1 ; 

Vf d (p,q,k) = TC p (p,q,k) (12) 
A(*? AfJ . B(p 2 ,iJ, 2 )-B(q 2 ,fi 2 ) / 

{i> +i)ip + q)v 
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, Mp 2 


» 2 )- 


A(q 2 ^ 2 ) 




2(p 2 - 


-q 2 ) 1 


A{p 2 




A(q 2 ^ 2 ) 



[{p 2 - q 2 h» - (P 1 -4)(p + q)u] 

p 2 + q 2 , . 

(p 2 -q 2 ) 2 + (M 2 (p 2 ) + M 2 (q 2 )) 2 1 ; 

with the new parameter d to satisfy the conditions (i) and (ii). Again k denotes the gluon 
momentum and p and q the quark momenta. The anomalous dimension 5 of the ghost 
propagator is 5 = —9 N c / (44 N c — 8 Nf), i.e. the corresponding value at one loop order for 
N c colors and Nf flavors. The Abelian part of the vertex is given by the Curtis-Pennington 
(CP) vertex T^ p (p,q,k). 

From a systematic point of view the newly introduced parameter d in the non-Abelian 
part of the vertex is completely arbitrary. Our numerical results, however, will indicate that 
values around the somewhat natural choice d = match best with lattice simulations, see 
below. 

For comparison we will also employ the much simpler vertex T u = y^ bel \y^ abel with W^ abel 
given by eq. (jTTj) and 

V: bd (p,q,k) = Z 2 (/i 2 ,A 2 ) 7 , (14) 

where we have taken the bare Abelian vertex, 7^, multiplied with an extra factor of Z 2 . In 
Landau gauge this construction also satisfies the conditions (i) and (ii), as will be shown 
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in the next two subsections. Furthermore in the numerical treatment we will additionally 
employ a vertex where the last, transverse term in Eq. (|13|) is left out, i.e. a generalized 
Ball-Chiu (BC) vertex. In Landau gauge such a vertex also satisfies the conditions (i) and 
(ii). 



Multiplicative renormalizability of the quark equation 

To proceed we first substitute the simpler vertex ansatz (f*HJ) into the quark equation (fTJ) . 
By taking the Dirac trace once with and once without multiplying the equation with p 1 we 
project out the mass function M(p 2 ) and the vector self energy A(p 2 ). To ease notation we 
will use in the following the abbreviations x = p 2 , y = q 2 and z = (p — q) 2 for the squared 
momenta, also s = fi 2 for the squared renormalization point and L = A 2 for the squared 
cutoff of the theory. We arrive at 



M(x)A(x,s) = Z 2 (s } L)m (L) + ^-^- / <Pq\ \' Z 2 (s, L) A (y, s) x 

37T 3 J \z(y + M 2 (y)) 



2d-d/8 

G(z,s)Z 3 (s,Ly 



3M(y) } , (15) 



(Z( Zl s)Z 3 (s,L)) d 

^ s) = Z 2 (s, L ) + ^1J Z 2 ( S , L) A-i(y, s) x 

/ \ -2d-d/8 

(G(z,s) Z 3 (s,L)j f +y ( X -y f 



(Z(z,s)Z 3 ( Sl L)) d V 2 2z 

(16) 

where we have used the following expression for the running coupling a in Landau gauge [^] 
a(x) = — Z(x, s)G 2 (x, s) = a(s) Z(x, s) G 2 (x, s) . (17) 

4:71 

From eqs. (|T5|) and (|T6|l it is clear that the choice of the parameter d = is special since 
then the only input from the Yang-Mills sector is the running coupling. 

The behavior of eqs. (|15jl. (|TH|) under renormalization can be explored by changing the 
renormalization point s = fi 2 to a new point t = v 2 . We first note that the factor stemming 
from the non-Abelian part of the quark-gluon vertex is not affected by such a change: 

2d-d/S / \ -2d-d/5 



(G(z,s)Z 3 (s,L)) (G(z,t)Z 3 (t,L) 



(Z(z,s)Z 3 (s,L)) d " (Z(z,t)Z 3 (t,L)) d 

This can be seen easily with the help of the relations 

G Q (x,L) = G(x,s)Z 3 (s,L), (18) 
Z Q (x,L) = Z(x,s)Z 3 (s,L), (19) 

between the unrenormalized and renormalized ghost and gluon dressing functions. Further- 
more, the running coupling at(z) is independent of the renormalization point, for a detailed 
discussion of this property using the expression (|T7|) see ref. 0. 
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A(t,s) 
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A(x,t) 
■ *- 



FIG. 2: Sketch of a finite renormalization from a perturbative point s to a non-perturbative point 
t for the vector self-energy A. 



From eq. © we infer 

Z 2 {t,L) A~\x,t) = Z 2 (s,L) A' 1 
With the renormalization condition A(t,t) = 1 we have 

Z 2 (t,L) = Z 2 (s,L) A~\t, s) 

and subsequently 

A(x,t) = A(x, s) A" 1 ^, s) . 



x,s). (20) 



(21) 



(22) 

Substituting eqs. (}2"T]) and into the Dyson-Schwinger equations (fTSJl we find the mass 
function M(x) to be independent of the renormalization point, i.e. condition (i) is satisfied. 
Note that without the extra factor of Z 2 in the Abelian part of the vertex (|14|) we would 
violate this condition. 

Before we examine the case of the more sophisticated Curtis-Pennington type vertex (|13|) 
two remarks are in order. First, according to perturbation theory we have A(x — > oo, s) — ► 1 
and Z 2 (s,L) — > 1 for large renormalization points s. However, this is just a special case of 
the general relation 

A(x -> oo,s) -> Z 2 (s,L) (23) 

which can be inferred from eqs. ()21j) and ()22j). In Fig. Elwe sketch the vector self energy 
renormalized at two different points s and t, with s in the perturbative and t in the non- 
perturbative region of momentum. Of course, this mechanism will be found again in the 
numerical results, see below. Second, the appearance of the ghost and gluon renormaliza- 
tion factors Z 3 (L) and Z 3 (L) in the interaction kernel of the quark equation is due to the 
non-abelian part of the employed quark-gluon vertex ansatz. Certainly, the renormalized 
functions M(x) and A(x, s) should not depend on the cutoff of the integral. The balance 
of cutoff dependent quantities in the equation is controlled by various factors of Z 2 (L) and 
Z m (L). We have to take care not to disturb this balance by the vertex ansatz. _ Thus the 
non- Abelian part of the quark-gluon vertex contains such powers of Z 3 (L) and Z 3 (L) that 
the cutoff dependence of these quantities cancel. This can be easily checked using the scaling 
behavior 



Z 3 (s,L) 



Z 3 (s,L) 



a{L) 



ays) 



a{L) 
a(s) 



(24) 



S 



of the renormalization factors for L — > oo and the relation 7 + 25 + 1 = for the anomalous 
dimensions of the gluon and the ghost, respectively. 

Along the same lines as for the bare vertex construction we prove condition (i) for the 
Curtis- Pennington type vertex. Plugging eqs. (fTTj). (fT3j) into the quark equation (JTJ) and 
projecting onto M(x) and A(x) we arrive at 



M(x)A(x,s) = Z 2 {s,L)m {L) 



a(z) 



G(z,s)Z 3 (s,L) 



-2d-d/8 



z(y + M 2 (y)) (Z{z,s)Z 3 (s,L)) c 



xA- l {y,s) -(A(x,s) + A(y,s))M(y) 

+i(AA M(y) - AB) f-z + 2{x + y) - 



+-(A(x, s) - A(y, s))M{y)Q{x, y){x - y) 



(25) 



A(x, s) = Z 2 (s,L) 



> dq 



a(z) 



G(z,s)Z 3 {s,L) 



-2d-d/5 



xz(y + M 2 (y)) (Z(z, s)Z 3 (s, L)) c 



xA ifas) 



x + y (x - y) 2 \ A(x,s) + A(y,s) 



/ AA 



2z 



(x + y) + AB M(y) 



Z - + (x + y)- ix ~ y)2 
2 v yj 2z 



+ l{A{x, s) - A(y, s))M(y)n(x, y) _ Z X _A 



Here we have used the abbreviations 



AA 
AB 

Q(x,y) 



A(x,s)-A(y,s) 

x — y 
B(x,s)-B(y,s) 



x-y 



x + y 



{x - y) 2 + (M 2 (x) + M 2 (y)) 2 ' 



(26) 



With the help of the relations (}2~T]) and (J2~2~j) we find again the quark equations (J23j) . (J2TI)) to 

be consistently renormalized. Employing 



B(x,t) = B{x,s) A~~ l (t, s) 
it can be seen directly that the mass function 

M(x) = B(x,s)/A(x,s) 



(27) 



(2* 



is independent of the renormalization point. (The same is true for a Ball-Chiu type vertex, 
which is the Curtis-Pennington construction (jlHjl without the transverse term proportional 
to Q(x,y). Note that in different gauges than Landau gauge only the Curtis-Pennington 
construction would satisfy condition (i), similar to QED |23||.) 
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Ultraviolet analysis of the quark equation 



In this subsection we will show that the ansdtze ()13|) and ()14j) for the quark-gluon vertex 
both lead to the correct perturbative limit of the quark mass function M(x). We first 
examine the case of the bare vertex construction, eq. (|15|) . 

The ghost and gluon dressing functions G and Z are slowly varying for large momenta ac- 
cording to their perturbative limit. For loop momenta y larger than the external momentum 
x we are therefore justified to employ the angular approximation G(z), Z(z) — > G(y), Z(y), 
see ref. Q. Furthermore, there is a region xq < y < x where the approximation 
G(z), Z{z) — * G(x), Z[x) is adequate. We are then able to carry out the angular integrals in 
eq. (fTo^) . If we additionally take the external momentum x to be large enough all masses in 
the denominators become negligible since the integral is dominated by loop momenta y ~ x. 
We then obtain 



M(x)A(x,s) 



Z 2 (s,L)m (L) 
Z 2 (s,L) a(x) 



+ 



+ 



7T 



X 



dyZ 2 (s,L)A (y,s) 



G(z,s)Z 3 (s,L) 



-2d-d/S 



(Z(z,s)Z 3 (s,L)) d 

-2d-d/8 



■M(y) 



Z 2 {s,L) f a(y) 



7T 



dy- 



Z 2 (s,L) A-\y,s) 



G(z,s)Z 3 (s,L) 



[Z(z,s)Z 3 (s,L)) c 



-M(y) , 
(29) 



where the integral from y = to y = xq has already been neglected. 

For large momenta y > xq the wave function renormalization A' 1 and the renormalization 
factor Z 2 cancel each other according to eq. ()23|) . The ultraviolet limit of the ghost and gluon 
dressing functions from their respective DSEs in quenched approximation has been discussed 
in ref. |4j and found to be in agreement with resumed perturbation theory to one-loop order. 
As will be seen in the discussion beyond quenched approximation below (see eqs. (JHHJ) ) , we 
also obtain the correct anomalous dimensions for the ghost and gluon dressing functions G 
and Z in the case of Nf ^ 0. We thus use the general perturbative limit 



G(z) 
Z(z) 



G(s) 
Z(s) 



^log (^-J +1 
^log(^) +1 



(30) 



with to = j3 a(s) / (An) = (lliVc — 2Nf)a(s)/ (12tc). If we additionally substitute the scaling 
behavior of the renormalization constants Z 3 and Z 3 , eqs. (|24|) . and exploit the relation 
7 + 25 + 1 = we arrive at 



M(x) = m (L) + 



1 a(x) 



ix x 



dyM{y) + 



TT 



(31) 



A) 



This well known equation describes the ultraviolet behavior of the quark mass function, 
see e.g. ref. (24j ]. Employing the perturbative form of the running coupling, 



QMS 



LO log ( - 

S 



1 



-1 



(32) 
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one obtains in the chiral limit, mo(L) = 0, the so-called regular asymptotic form 



M(x) = ^ 



3 x(iln(x/A^ D )) 



1—7" 



(33) 



Here (^^I/) denotes the chiral condensate which is discussed in more detail in the next 
subsection. In the case of non- vanishing bare quark mass, mo(L) ^ 0, the equation (j31j) is 
solved by the irregular asymptotic form, 



M{x) = M{s) 
In this case we furthermore find 



(^) + I 



12 



-"Ira 



llN r - 2N f ' 



m (L) = M(s) 



f 

wlog ( ^ ) I 



"7m 



(34) 

(35) 
(36) 



in accordance with perturbation theory. 

We thus have shown that the bare vertex construction (J14|) admits a solution for the 
mass function M(x) which has the correct perturbative behavior for large momenta. A 
similar analysis is possible for the DSE with the Curtis-Pennington type vertex, eq. (|25|). 
As the vector self energy goes to a constant in the limit of large momenta, eq. ()23|) . all 
terms proportional to A(x) — A(y) are suppressed in this limit. Furthermore, according to 
the perturbative expression ()34j) the A5-term contributes at most subleading logarithmic 
corrections in eq. ([25)1 . The first term in the brackets reduces to the bare vertex form because 
A(x, s) ~ A(y, s) for large momenta x, y. Thus we obtain the same ultraviolet limit from 
eq. (|25|) as for the bare vertex construction. This is certainly also the case if a Ball-Chiu 
type vertex is employed. 



THE QUARK PROPAGATOR IN QUENCHED QCD 

In this section we will compare quark propagator results in quenched approximation for 
three different vertex types, which share the non-Abelian part proposed in eq. (jllj) but differ 
in their Abelian parts. We will employ the bare vertex, eqs. (|14|). and the Curtis-Pennington 
(CP) type vertex, eqs. (|13|) . Furthermore we use a Ball-Chiu (BC) type construction, which 
employs only the first three terms of the CP-vertex. In Landau gauge all these vertex ansatze 
satisfy the conditions (i) and (ii) formulated above eq. (Jll)) . In order to compare the different 
vertex types on a quantitative level we will calculate the pion decay constant f n and the 
chiral condensate from the respective solutions for the quark mass function. 

Before doing so we have to specify the effective quark interaction as input from the Yang- 
Mills sector. 



Effective quark interaction 

Via the ansatz of the quark-gluon vertex the effective quark interaction depends on the 
quark propagator functions A(x, s) and M(x) themselves. These functions will be determined 



11 





FIG. 3: Shown are the results for the gluon dressing function Z, the ghost dressing function G and 
the running coupling a. The two sets of fit functions are given in eqs. (|37|) and (|38|). The /^-function 
corresponding to our DSE-solution is compared to the one- and two-loop expressions as well as to a 
polynomial in a. 



self-consistently in the process of the solution the quark DSE. As further input the running 
coupling and, if d ^ 0, the gluon and the ghost propagator are needed. These will be taken 
from the results of ref . 0] ■ 

The DSEs of the Yang-Mills sector are scale-independent and the underlying scale will 
be generated through dimensional transmutation during renormalization. To translate this 
scale into physical units the scale Aqcd is, in quenched approximation, determined from the 
large-momentum behavior of the running coupling, or more precisely, by requiring a(M§) = 
a((91.187 GeV) 2 ) = 0.118. Technically this is achieved by fitting the functional form of the 
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running coupling. We employ two different fit functions [25] for the running coupling a(x): 

Fit A: a(x) = — r -, — t-t^ — — —, — -r—^ rr-r , (37) 

ln[e + ai (x/A 2 QCD )^ + hix/Al^} 

Fit B: a(x) = : — —-7, rr a a(0) + 

a + (x/A 2 QCD f [ 



(3* 



The value a(0) = 8.915/iV c is known from the infrared analysis. In both fits the ultraviolet 
behavior of the solution fixes the scale, Aqcd = 0.71 GeV. Especially the fit B is very robust 
leading to a change in Aqcd of less than one percent when different fitting regions are 
chosen. Note that we have employed a MOM scheme, and thus Aq C d has to be interpreted 
as A^/qm, i.e. this scale has the expected magnitude [52j] . Fit A employs the four additional 
parameters: a\ = 1.106, 02 = 2.324, b\ = 0.004, 62 = 3.169. Fit B has only two free 
parameters: a = 1.020, b = 1.052. 

In Landau gauge the gluon and ghost propagators are defined by 



D M = (V - P -f) '-f- , (39) 
G(p 2 ) 

D G (p) = ^ . (40) 

The gluon and ghost dressing functions, Z{p 2 ) and G(p 2 ), respectively, are then described 
by 



R(x) 



e(x/A 2 QCD r + d(x/A 2 QCD r 
l + c(x/A 2 y + d{x/A 2 ) 2 - 



V QCD) ~ u W 1Y QCD 
Z(.r) = [^1) 1+2S R 2 ( X ), 

g{x) = (441 fl-'c*). ( 4i ) 



a(/i) 

where x = p 2 and c, d are fitting parameters for the auxiliary function R(x). They are given 
by c = 1.269 and d = 2.105. Recall that the anomalous dimension 7 of the gluon is related 
to the anomalous dimension 5 of the ghost by 7 = —1 — 26 and 5 = —9/44 for the number 
of flavors Nf = 0. 

Whereas Fit A is better in the region 0.3 GeV 2 < x < 1 GeV 2 where a is strongly rising, 
Fit B is slightly better in the region 1 GeV 2 < x < 10 GeV 2 . As can be seen in Fig. Elboth 
fits works very well and will be used as input below. 



Pion decay constant, chiral condensate and quark masses 

A correct calculation of the pion decay constants involves the pion Bethe-Salpeter am- 
plitudes including the subleading components, see e.g. refs. jU!^]. Apart from the dressed 
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quark propagator the Bethe-Salpeter equation involves couplings between quarks and gluons. 
On the level of the quark DSE we have substituted the full quark-gluon vertex by an vertex 
ansatz. However, at present it is only known for certain cases how such a vertex ansatz in 
the quark DSE translates to the corresponding quark-gluon coupling in the Bethe-Salpeter 
equation |2^|. No method is known up to know to derive the corresponding Bethe-Salpeter 
equation for dressed quark-gluon vertices as the BC- or CP-vertex constructions. 
We thus have to rely on the approximation |3(| 

f 2 N n f , M(y)A- 1 (y,8) ( y dM{y) \ 

which incorporates only the effects of the leading pion Bethe-Salpeter amplitude in the chiral 
limit |3llj. From a comparison of the relative size of the amplitudes in model calculations 
(HI l33l] one concludes that the approximation (J4*2~j) should lead to an underestimation of /„. 
by 10-20 %. 

The renormalization point independent chiral condensate, (\£ / \I / ), can be extracted from 
the ultraviolet behavior of the quark mass function in the chiral limit (c.f. eq. (|33|)): 



27T 2 7m -(jng) 

x (lln{x/A 2 QCD )) 



m{x) ^ *\' m - d :*^ (43) 



Hereby Aqcd is to be taken from a fit to the running coupling, c.f. eqs. fl3Tj) . (|3*8J). 

The renormalization point dependent chiral condensate (^ rx I / ) M can be calculated via |2 

J^S ch (q 2 } s), (44) 

where the trace is over Dirac indices, S c h is the quark propagator in the chiral limit and the 
squared renormalization point is denoted by s = fi 2 . To one- loop order both expressions for 
the condensate are connected by 

(**>m= Qm^/Aqcd)) 7 " 1 (45) 

with 7 m being the anomalous dimension of the quark mass function. 

For the calculation of the chiral condensate we first have to determine the mass renor- 
malization constant Z m (s,L). Recall the formal structure of the mass equation (|15jl . which 
is given as 

M(x)A(x, s) = Z 2 (s, L) Z m (s, L) m R {s) + Z 2 (s, L) U M (x, s) , (46) 

where n M (x, s) represents the dressing loop. In order to extract Z m (s, L) from this equation 
we have to clarify the meaning of m^(s) which is related to the unrenormalized mass by 

m {L) = Z m (s, L) m R (s), (47) 

c.f. eq. (0). Evaluating eq. ()46j) at the perturbative momentum x = s the matter seems 
clear. We achieve consistency with eqs. and (JHHJ), if 

m R (s) = M(s) , (48) 



Z m (s, L) 



-■log ( I ) 4 1 



"7m 



(49) 



14 



which is indeed the correct perturbative scaling of the renormalization constant Z r , 

Certainly one could implicitly define the finite parts of Z m such that the relation (}4*%)) holds 
in general for all renormalization points s. Then the parameter m R in the renormalized QCD- 
Lagrangian would already know about dynamical symmetry breaking. However, as the mass 
parameters of QCD are supposed to be generated in the electroweak sector of the standard 
model one could equally well argue that it is more systematic to exclude the effect of mass 
generation by strong interaction from m R . 

In our numerical calculations we will choose s to be sufficiently large, therefore eq. 
is valid anyway. Then Z m is determined by 



Z m (s, L) 



M(x) A(x, s) - Z 2 {s, L) Il M (x,s) 

(50) 



Z 2 (s,L)M(s) 
U M (s,s) 



Z 2 (s,L) M(s) 



For the last equation we have set x = s and have used the renormalization condition A(s, s) = 
1. 

In the numerical calculations we have to specify the masses m R (s) as input. Choosing a 
perturbative renormalization point s allows one to evolve the masses m R (s) to a different 
scale t by 



( H S / A QCD 
\Mt/A$lCD 



7m 



m R (t) = m R (s) - / /A r" • (51) 



For t = (2 GeV) 2 typical values for the masses of the light quarks are given by the Particle 
Data Group Q: 

^(m u + m d ){2 GeV) w 4.5MeV, m s {2 GeV) w lOOMeV. (52) 
We will use similar masses in our calculations. 



Renormalization scheme and numerical method 

In the quark equation we employ a MOM regularization scheme similar to the one used for 
the ghost and gluon equations in refs. 0,0]. For a fermion DSE this technique has already 
been used in quenched QED 4 [35]. The formal structure of the quark equation is given by 

A(x, s) = Z 2 {s, L) + Z 2 (s, L) U A {x, s) , (53) 
M{x)A{x, s) = Z 2 (s, L) Z m {s, L) m R {s) + Z 2 {s, L) n M (x, s) . (54) 

We eliminate Z 2 from the first equation by isolating it on the left hand side and subtracting 
the same equation for x — s. With 

' 1 + TT~ r U A (x,s), (55) 



Z 2 (s,L) A(x,s) A(x,s) 

we then have 

1 1 



A(x,s) A(x, s] 



U A (x,s) + U A (s,s), (56) 
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using the renormalization condition A(s, s) = 1. In each iteration step we determine A(x) 
from eq. (JSTIj) and subsequently Z 2 from eq. (|55|). As a numerical check we determine Z 2 at 
different momenta x = p 2 . In our calculations we find Z 2 to be independent of p 2 to a high 
accuracy. For the mass function M(x) we use 

M(x)A(x, s) = Z 2 (s, L) Tl M (x, s) (57) 

in the chiral limit and the subtracted equation 

M{x)A{x, s) = M{s) + Z 2 {s, L) U M {x, s) - Z 2 {s, L) U M {s, s) (58) 

if chiral symmetry is broken explicitly, i.e. mo 7^ 0. 

For the numerical iteration we employ a Newton method and represent the dressing 
functions A(x) and M(x) with the help of Chebychev polynomials. Furthermore, we use 
a numerical infrared cutoff e, which is taken small enough for the numerical results to be 
independent of e. Numerical difficulties arise in the case of the Curtis-Pennington type 
vertex and even more for the Ball-Chiu construction. If the external momentum x and the 
loop momentum y are both small and close to each other then the derivative-like terms 

AA = Ajx, s)-A(y,s) AB = B(x,s)-B(y,s) 

x — y x — y 

are hard to evaluate accurately. Although the functions A(x) and B(x) are constant in the 
infrared and consequently should have derivatives close to zero one encounters large values 
for A A and AB due to numerical inaccuracies in A and B. In order to evaluate A A and 
AB much more precisely at small momenta we fit the expressions 

M^) = - y°>*l . , B(x,s) = - Vh , (60) 

with the parameters a 1; a 2 , &i and b 2 to the numerically evaluated functions. The scale 
Aqcd = 0.71 GeV has been determined from the fits to the running coupling already in the 
last subsection. 

For x — y smaller than a suitable matching point we calculate the terms AA and AB 
from the fits. This procedure eliminates the numerical errors in the derivative terms and 
smoothes the numerical results considerably. In the case of the Ball-Chiu type vertex the 
iteration process does not converge unless we use these fits. 

The renormalization condition employed in the ghost-gluon system of equations is 
G 2 (s)Z(s) = 1 with a(s) = 0.118 at the squared renormalization point s = a 2 . Furthermore 
we choose a transversal tensor to contract the gluon equation, c.f. ref. The physical 
scale in the quenched calculations is taken directly from the Yang-Mills results, i.e. we use 
the experimental value a(M§) = 0.118 of the running coupling at the mass of the Z-boson 
to fix the scale. 



Numerical results 

In Fig. E] we give our numerical solutions for the quark mass function and the inverse 
vector self energy in the chiral limit, employing 'Fit A' for the effective quark interaction. 
We compare results obtained with five different ansdtze for the quark-gluon vertex. For the 
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FIG. 4: The mass function M(x) and the inverse vector self energy 1/A(x) of a chiral quark are 
shown. We compare the results for five different vertices with lattice data taken from refs. [l^. 





M(0) 


U 


(-(§^))V3 


(-(**))V3 


A-\0,M 2 z ) 


a\ 


a.2 


bi 


b 2 




[MeV] 


[MeV] 


[MeV] (calc.) 


[MeV] (fit) 












bare vertex 


177 


38.5 


162 


160 


0.733 


3.05 


0.99 


0.06 


1.00 


CP d=5 


150 


50.5 


223 


225 


0.910 










BC- vertex 


293 


62.6 


276 


284 


0.523 


1.10 


0.99 


0.29 


0.92 


CP-vertex 


369 


78.7 


303 


300 


0.634 


0.83 


0.99 


0.20 


1.00 


CP d=0.1 


464 


87.5 


334 


330 


0.501 


0.79 


0.99 


0.34 


0.95 



TABLE I: The mass M(0), the pion decay constant f n calculated with eq. (|42|). the renormalization 
point independent chiral condensate calculated with eqs. (|44j) and (|45j) . and the condensate obtained 
by fitting the expression ()43j) to the chiral mass function in the ultraviolet for all five vertex types. 
Recall 5 = —9/44 in quenched approximation. If not stated otherwise the parameter d in the vertex 
construction is taken to be d = 0. For the case of the CP-vertex with d = 5 we did not get good 
fits in the infrared. 



generalized CP- vertex we investigate the 'natural' case d = 0, the value d = 5 = —9/44, 



already adopted in refs. |20j, |21|, and the value d = 0.1. Furthermore we employed the 
bare vertex construction and a Ball-Chiu type vertex. The corresponding masses at the 
momentum p 2 = 0, the pion decay constant the renormalization point independent chiral 
condensate and the fit parameters for the functions (J6U)) are displayed in table 

The numerical results for the mass function all have a characteristic plateau in the infrared 
and show the regular asymptotic behavior for large momenta, c.f. eq. (}4"3*|) . The bare vertex 
construction and the CP type vertex with d = 5 both generate masses much smaller than 
typical phenomenological values of 300 — 400 MeV. The BC- and the CP-type construction 
with d = provide good results, whereas the choice d — 0.1 leads to a somewhat large mass. 
The lattice calculations taken from refs. 0] (overlap fermions) and 0| (improved staggered 
action, Asqtad) favor masses around 300 MeV with the caveat that they are obtained by an 
extrapolation from sizeable bare quark masses to the chiral limit. The numerical solutions for 
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FIG. 5: Results for three different forms of the running coupling in the quark equation: The running 
coupling calculated in ref. 0] and the two fits given in eqs. (J37|). J3HJ). 



the wave function renormalization 1/A can be seen in the right diagram of Fig. 0] Whereas 
the ultraviolet asymptotic behavior of all vertex constructions is similar we observe sizeable 
differences for small momenta. Again the bare vertex construction and the CP-vertex with 
d = 5 are clearly disfavored by the lattice data. 

Our approximate calculation of the pion decay constant should underestimate the exper- 
imental value f n = 93MeV by 10-20 %, c.f. the discussion below eq. We thus have 
best results for the CP- vertex construction with d = and d = 0.1. Furthermore we obtain 
very good agreement between the two different methods to extract the chiral condensate. 
Compared to the phenomenological value (— (\E , \E')) 1 / 3 « 250 MeV most of our results are 
larger. 

Apart from the case d = 5 we obtain very good fits for the scalar and vector self energy, 
A(x) and B(x), for small momenta. The results for the parameters 01,02,61,62 in the 
fit functions given by eq. (jfjtjj) can be found in table |TJ It is interesting to note that the 
exponents 02 and 62 are found to be very close to one. Such a behavior could indicate a 
simple underlying functional form of the quark propagator. This will be explored in future 
work by a numerical continuation of our results to negative p 2 , i.e. timelike momenta. 

Fig. El compares results for the bare vertex and the CP-type construction with d = for 
three different forms of the running coupling in the interaction kernel of the quark equation. 
The two fit-functions, 'Fit A' and 'Fit B', have been given in eqs. (137)) . ()38j) . Furthermore 
we used the running coupling calculated from the quenched ghost and gluon DSEs in ref. . 
Although there is the (presumably) artificial bump at p 2 ~ 0.1 GeV 2 in the running coupling, 
the mass functions obtained from the DSE-result and from 'Fit A' are virtually indistinguish- 
able. 'Fit B', however leads to somewhat smaller masses. This observation suggests that 
nearly all the dynamically generated mass is produced from the integration strength above 
p = 500 MeV, indicated by the vertical line in the plot of the running coupling. This is a 
favorable result as it would have been very unsatisfying if the artificial bump contributed a 
considerable amount to the quark mass function. 
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FIG. 6: These diagrams show our results when three different bare quark masses are employed. In the 
diagram on the right small quark masses correspond to small values for X/A in the infrared. 



Finally we observe the effects of explicit chiral symmetry breaking in the plots of Fig. |UJ 
We give results for three different quark masses, m(2 GeV) = 5 MeV, m(2 GeV) = 100 
MeV and m(l GeV) = 1000 MeV. These values correspond roughly to the ones given by 
the Particle Data Group for the up/down-quark, the strange-quark and the charm-quark 
|34| . For small momenta we note again that the dressed vertex generates more mass in the 
quark equation than the bare vertex construction. This effect becomes much less dominant 
for the heavy quarks, where more and more of the infrared mass stems from explicit chiral 
symmetry breaking and not from dynamical mass generation. Furthermore in accordance 
with the analytical determination we observe the same ultraviolet behavior of the mass 
function for both vertex constructions. 



THE QUARK-LOOP IN THE GLUON DSE 

In this section we focus on the inclusion of the back-reaction of the quarks on the ghost- 
gluon system, i.e. we will finally solve the complete set of mutually coupled Dyson-Schwinger 
equations for the quark, gluon and ghost propagator. To this end we incorporate the quark- 
loop in the truncation scheme for the gluon DSE which has been developed in ref. pj. 

The formal structure of the gluon equation is given by (c.f. fig. 

[d{p%1 = z 3 [dW(p)]-J; + nf u ° st ( P ) + ngr(p) + n*r fc (p). (6i) 

The contributions from the ghost- and gluon- loop, H 9 Jl° st (p) and U 9 ^ on (p), are treated in 
detail in ref. j^|. The contribution of the quark- loop to the gluon equation is given by 

n*T*(p) = Zif I d'q Tr { 7m S(q) T v {q,p) S(k)} , (62) 
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FIG. 7: Diagrammatical representation of the quark-loop in the Dyson-Sch winger equation for the 
gluon propagator. 

where p is the external gluon momentum and q and k — (q — p) are the momenta of the two 
quarks running in the loop, (cf. fig. EJ). The trace is over Dirac indices. 

In eqs. fjl lH13j) we have proposed an effective quark-gluon vertex T u with Abelian and 
non-Abelian parts such that the quark equation is multiplicatively renormalizable and one- 
loop perturbation theory is recovered for large momenta. However, this construction is 
not capable to account as well for the one-loop behavior of the unquenched gluon equation 
unless we switch the momentum arguments of the non-Abelian part w^ abe \ eq. fTTf • from 
the gluon momentum to quark momenta. In the quark equation such a change of momentum 
arguments would either break Lorentz symmetry by preferring one quark line of the quark- 
gluon vertex or changes the ultraviolet behavior of the quark equation. We therefore have to 
use different momentum assignments for the quark-loop and the quark equation. Certainly, 
this is a deficiency which has to be resolved by a more elaborate vertex construction in future 
work. The aim of the present study, however, is to present an effective construction which 
captures essential properties of the theory. 

Taking care of symmetries we propose the following ansatz for the non-Abelian part of 
the quark-gluon vertex in the quark-loop: 



/ ~ \ -d-d/(2S) / _ \ -d-d/(28) 

- (G(y)Zs(L)) (G(z)Z 3 (L)) 
q P (Z(y)Z 3 (L)f 2 (Z(z)Z 3 (L)f 2 



(63) 



Here x = p 2 is the squared gluon momentum, y = q 2 and z = k 2 = (q — p) 2 are the squared 
quark momenta, and L = A 2 is the squared cutoff. The Abelian part V^ bel of the vertex, 
given in eqs. ()13I14|) . is symmetric with respect to the quark momenta as well. 

Plugging the Curtis-Pennington type vertex into the quark-loop and contracting the free 
Lorenz-indices with the tensor (c.f. the treatment of the gluon DSE in ref. 

P#(P) = ^-C^, (64) 
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we obtain 

TT 9 2 N fr7 f j4 G(y) Gjz) (G(y) G(z) Zg(L))-^/^ 

IWfc - (2?r)4 ^2 y fl ^ + M 2 (y) z + M 2( z ) Z(z ) 2|(L))d/2 

xA- 2 (y)A- 2 (^) ( ^(y) + ^ y , z)A(y)A(z) + W 2 (s, y, z)£(y)B(*)) 

(W 3 (a;, y, ^(y)^) + W 4 (x, y, z)B(y)B(z)) 
(W 5 (x, y, z)A(y)B(z) + W 6 (z, y, z)B(y)A(z)) 
(W 7 {x, y, z)A(y)A(z) + W 8 (x, y, z)B(y)B(z)) 



\ 2 




, A(y) - 


A(z) 


2(3/- 


-z) 


, B(y) - 


B{z) 



y - z 

(A(y) - A(z))(y + z) 
2((y - z) 2 + (M 2 (y) + M 2 (z)f 



(65) 



with the kernels 



Trr / , (z 2 /2-C 2Cy\ 2 (2-0y CV 2 ^ 

W a (x,y,«) = ^ZO, (67) 
3a; 

w( \ C ^ 2/l + C^CyV A, (2C-6)y Cl/ 2 



3x 2 \ 3x 3x 2 J \3 3a; 3a; 2 

y _ (C + i)y 2 + Cy 3 

3 3x 3x 



+^-^r^ + S' ( 68 ) 



u/r ^ ~2C^ 2 , ( 4 4(y\ 2 4y 2(y 2 .... 

3x z \3x 3x z J 3 3x 3x z 



= 3^-H^ + ^ + 3 + ^^ + 3^' (70) 



M/r ^ ^ 2 , 2(y^ , 1 , (-C-l)y , Cz/ 2 , 7n 

2 2 

lVr(x,y,2) = - Z - + z + V -- yi (72) 

W 8 (ar,y,z) = 2 (-- + ^) . (73) 
V x xJ 

Note that the symmetry factor 1/2 and a factor l/(3x) from the left hand side of the gluon 
equation have been absorbed in the kernels. From this expression the corresponding one for 
the bare vertex construction can be read off easily by setting = and replacing the 

remaining factor (A(y) +A(z))/2 in eq. by unity. 



Ultraviolet analysis of the quark-loop 

It is long known that the introduction of a cutoff A in the gluon DSE results in artificial 
quadratic divergencies due to the violation of gauge invariance. Certainly, to recover the 
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correct perturbative limit of the gluon propagator such terms have to be removed from the 
gluon equation by a suitable regularization procedure. Quadratic divergencies only occur in 
the part of the inverse gluon propagator proportional to 8^. Therefore one way to eliminate 
the quadratic divergencies is to project onto the part proportional to p^p u |36[ by choosing 
( = 4 in the projection tensor (|64j) . 

Another unambiguous way is to subtract the quadratically divergent terms from the 
kernel by hand. This procedure is valid for general ( and allows to estimate the influence of 
spurious longitudinal terms in the right hand side of the gluon equation on the solutions by 
varying the parameter (. Certainly in a perfect truncation scheme the solution should be 
independent of (. 

In ref. 0] it has been described in detail how to remove the quadratic divergencies in 
the ghost- and gluon-loop. Therefore in the following we concentrate on the quark-loop. 
To identify the divergent terms we expand the dressing functions in the integrand of the 
quark-loop around large loop momenta y with the difference (z — y) still larger than any 
quark mass. 

To leading order this expansion amounts in the replacements 

G(z) - G(y), 
A(z) -> A(y), 



A{y) - A{z) 
y- z 

B(y) - B(z) 
y-z 

(A(y) - A(z))(y + z) A'{y){y + z 



B'(y), 

(74) 



2((y - z)* + (M*(y) + M^z))*) 2(y - z) ' 

with the derivatives A' and B' . Note that the first two equations are identical to the angular 
approximation employed previously in the ultraviolet analysis of the quark equation. For 
large momenta x and z the denominators in eq. (|fi5|) simplify and the angular integrals are 
trivially performed using the integrals given in appendix A. We arrive at 



KL* = / d V G\y) *-'(„) x 



Kyj \ 3y 3x 3xy yy! 

A>(y) fl + -2(4-Q g+ M + 2(4-0 W 
^3 3x \3y 3x J 

2 4-C\ . A'{y) fly , , 4„ 2 , 



+B\ y )M(y) y- - + ^ ^ - 1 + -tt'W j | . (75) 

Keeping in mind a factor (1/y) hiding in the derivatives we identify three quadratically 
divergent terms: (4 — ()/3x in the second line, —2(4 — ()y/3x in the third line and 4y/x in 
the last line. The first two of them are artefacts of the regularization and will be subtracted 
from the kernels. However, we encounter the additional (^-independent quadratic divergent 
term 4y/x originating from the transverse part of the Curtis- Pennington vertex. Such a term 
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is already known from corresponding studies in QED [37[. Although first suggestions have 
been made how the Curtis- Pennington vertex should be modified to avoid this problem 0], a 
convincing solution has not been found yet. In our study we therefore choose the pragmatic 
strategy of subtracting this term by hand together with the other quadratically divergent 
parts. 

Moreover we subtract all further terms proportional to (4 — (). Although these terms 
are not quadratically divergent they are artefacts of the regularization. We then obtain 
a (^-independent expression for the quark loop at large momenta. In ref. similar £- 
independent expressions for the ghost- and gluon-loop have been derived. We therefore 
arrive at a transversal right hand side of the gluon equation for large momenta as required 
in Landau gauge. 

Collecting all modifications together we have the new kernels 



Wi(x,y,z 
W 2 (x,y } z 
W 3 (x,y,z 

W 4 (x,y,z 

W 5 (x,y,z 

W 6 (x,y } z 

W 7 (x,y,z 
W 8 (x,y,z 



Wi(x,y,z) 
0. 

W 3 (x,y,z 

W 4 (x,y,z 

W 5 (x,y,z 

W 6 (x,y,z 

W 7 (x, y, z 
W 8 (x,y,z 



+ 



(y + z)(4-Q 

6x 



2zy{A - C) 
3x 

(y + z)(4-() 





3x 


(y + 


*)(4-C) 




6x 


(y + 


*)(4-C) 




6x 


(y- 


z)(y + z) 


X 



(76) 
(77) 
(78) 

(79) 

(80) 

(81) 

(82) 
(83) 



Note that the subtracted terms are chosen to preserve the symmetry of the kernels with 
respect to the squared quark momenta y and z. 

Without quadratic divergences we are in a position to extract the leading logarithmic 
divergence of the quark-loop. With modified kernels the ultraviolet limit of the quark-loop 
is given by 



jtUV 

quark 



9 2 N f 
16tt 2 ' 



n< v \-2d-d/6 



Z(y) d 
x -2 A'(y) 
3y 2 



\ + ^M\y) 
3 3y 



+B r (y)M(y)^-+^l 



-M\y) 

x 



34) 



Similar to the situation in the DSE for the quark mass function the leading ultraviolet term 
is the first term in the curly brackets. For the ghost and gluon dressing functions we employ 
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the perturbative ansatz 



G{z) = G(s) 
Z(z) = Z(s) 



UJ lot 



o;log(-) +1 



J5) 



and determine the anomalous dimensions 5 and 7 as well as the coefficient uj selfconsistently 
as follows. Substituting the ultraviolet limit of the vector self energy, eq. (j22J), and choosing 
the perturbative renormalization condition G(s) = Z(s) = 1 we arrive at 



u uv 



quark 



(P) 



2N t 



3(25 + IV 16vr 2 



-'loR ( — 

S 



2(5+1 



uj log ( - 

s 



2(5+1 



i6) 



Combining this expression with the results for the ghost- and gluon-loop from ref. |4| we 
obtain as ultraviolet limit of the gluon equation 



X\ 1 — T 

o;log( 7 j+l_ 



Z 3 + 



A^ 2 



+ 



27V^ 2 



96vr 2 cj(2<5 + 1) 48tt 2 u;(25 + 1) 48vr 2 cj(2<5 + 1) 

2(5+1 



-■log ( ~ ) + 



o;log(-) +1 



2(5+1 



57) 



The corresponding expression for the ghost equation reads 



x 



wlog(-) + 1 



64tt 2 cj(7 + 5 + 1) 



cj log 



L 



x 



7+5+1 



a; log 



x 



7+5+1 



The renormalization constants Z 3 (s,L) and Z 3 (s,L) cancel the cutoff dependence, i.e. the 
respective first terms in the curly brackets. Thus, the power and the prefactor of the second 
term have to match with the left hand side of the equations. This leads to the anomalous 
dimensions 



7 
6 



-13N C + 4N f 
22N C - AN f 

-9N C 
44A+ - 8iV 



59) 



/ 



which are related by 7 + 25 + 1 =0 and in accordance with one-loop perturbation theory 
for arbitrary numbers of colors N c and flavors Nf. For the coefficient uj one obtains 



uj = (11N C - 2iV/)a(s)/(127r) = #,a(s)/(47r). 



(90) 



When combined according to eq. (fTTjl our ghost and gluon dressing functions lead to the 
correct one- loop running of the coupling a(x) at large momenta. 
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Infrared analysis of the quark-loop 



The infrared analysis of the ghost anclgluon DSEs in a truncation with bare ghost-gluon 
vertex has been performed in refs. 0, 0, M To leading order the power law ansatz 

Z(x) ~ x 2K , G{x) ~ x~ K , (91) 

for the ghost and gluon dressing functions at small momenta x = p 2 has been employed. For 
a transverse projection tensor ie. £ = 1> one obtains k = (93 — VT201)/98 ~ 0.5954 
for the exponent of the dressing functions and subsequently the fixed point a(0) ~ 8.915/iV c 
for the running coupling in the infrared. 

These results have been obtained in a truncation where the ghost-loop dominates the 
gluon-loop in the infrared. Therefore in order to investigate the effects of dynamical quarks 
in the system we compare the infrared behavior of the quark-loop with the one of the ghost- 
loop. Substituting the ansatz (pTTj) into the gluon equation and calculating the ghost-loop 
along the lines of the infrared analysis given in refs. 0, 0| one finds the ghost-loop to be 
proportional to x~ 2k . For the quark-loop including the effects of dynamically generated 
quark masses we obtain 

KuarM-^ +2+ ^ & . (92) 

Therefore the quark-loop is suppressed for small momenta provided the parameter d fulfills 
the condition 

d < —. (93) 

K 

As we have k »s 0.5954 and 5 = — 1/4 for N c = 3 and Nf = 3 we find the condition d < 0.84, 
which is satisfied for all quark-gluon vertices employed in our calculation. From a numerical 
point of view we encounter serious instabilities in the quark equation once d is taken to be 
larger than d ~ 0.2. 

We conclude that with dynamically generated quark masses the quark-loop does not 
change the infrared behavior of the ghost and gluon dressing functions found in refs. 0, H| . 
In pure Yang-Mills theory as well as in QCD we thus have an infrared finite or vanishing 
gluon propagator and a ghost propagator which is more divergent than a simple pole. The 
Kugo-Ojima confinement criterion |38l.l39l| and Zwanziger's horizon condition [a El are both 
fulfilled not only in pure Yang-Mills theory (see e.g. refs. [25]) but also in QCD. 

The reasoning above reveals a selfconsistent picture valid for a small number of flavours: 
As has been demonstrated in our quenched calculations the combined strength of the ghost 
and gluon propagator generates sizeable dynamical quark masses in the quark DSE. In the 
unquenched case these masses suppress the quark loop in the infrared such that the ghost and 
gluon propagators are hardly changed in the infrared and in turn nearly the same amount of 
mass is generated in the quark equation as in the quenched case. This scenario is verified by 
our numerical calculations presented in the next section. On the other hand for a sufficiently 
large number of light flavours we expect a different selfconsistent picture to apply: Chiral 
symmetry should be restored. With vanishing quark masses the quark loop will contribute 
to the gluon DSE at small momenta and the infrared behaviour of the Yang-Mills sector will 
be changed. Especially it is expected that the value of the fixed point of the running coupling 
is decreased dramatically. This in turn drives the quark equation to the chirally symmetric 
solution with M(p 2 ) = 0. This second selfconsistent scenario as well as the situation under 
the presence of a small amount of explicit chiral symetry breaking are subject to future 
investigations. 
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M(0) 




h 


(_(^))l/3 
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(M z ) 


a MOM 
1V QCD 




[MeV] 


[MeV] 


[MeV] 




- 


[MeV] 


vertex 


qu. 


unqu. 


qu. 


unqu. 


qu. 


unqu. 


qu. 


unqu. 


qu. 


unqu. 


bare d=0 


177 


176 


38.5 


38.4 


160 


170 


0.118 


0.146 


710 


748 


CP d=5 


150 


133 


50.5 


46.3 


225 


230 


0.118 


0.146 


710 


746 


CP d=0 


369 


360 


78.7 


78.7 


300 


310 


0.118 


0.143 


710 


672 


CP d=0.1 


464 


437 


87.5 


85.5 


330 


330 


0.118 


0.142 


710 


640 



TABLE II: A comparison between the quenched (qu.) and unquenched (unqu.) results for the 
quark mass M(0), the pion decay constant /„-, the renormalization point independent chiral condensate, 
the running coupling at the mass of the Z-boson and ^qq^ for different vertices and values of the 
parameter d. The unquenched calculations are done for Nf = 3 chiral quarks. Furthermore we have 
8 = -9N C /(A4N C - 8N f ) = -0.25 in the present case. 



NUMERICAL RESULTS 

The numerical treatment of the integrals in the ghost and gluon equations has been 
described in detail in ref. |4(. The iteration process is done for the ghost-gluon system and 
the quark equations separately: we first iterate the Nf mutually uncoupled quark systems 
until convergence is achieved, feed the output into the ghost and gluon system, iterate until 
the ghost-gluon system converges, feed the output back into the quark equations and so on, 
until complete convergence of all equations is achieved. We renormalized at the point s = fi 2 
given by a(s) = 0.2 and used a transverse tensor to contract the gluon equation, ( = 1. 

In contrast to the quenched calculation we fix the physical scale of the system not by 
the condition a(Mf) = 0.118 but by adjusting the pion decay constant. This choice has 
an important advantage: Whereas the behavior of the running coupling for large momenta 
depends strongly on Nf the pion decay constant turns out to be almost independent of the 
number of flavors. In our quenched calculation we obtained f n = 78.7 MeV for the CP-type 
vertex with d = 0. Recalling that our approximation (}4"2*|) should lead to an underestimation 
of f n by about 10-20 % this value is in good accordance with experiment as f% xp ' = 93 
MeV. We therefore chose the scale in the unquenched calculations to lead to the same decay 
constant for the CP-type vertex with d — 0. 

In table |H] we compare results for the quenched and unquenched system of equations. 
The quark mass, the pion decay constant and the chiral condensate differ only slightly for 
each vertex construction respectively. The only sizeable difference occurs in the running 
coupling. As expected from perturbation theory the unquenched running for Nf = 3 results 
in larger values of the running coupling at p 2 = (Mz) 2 compared to the quenched case 
Nf = 0. We obtain a(M z ) ~ 0.140, which is somewhat larger than usually quoted values 
from experiment. However, such large values are not yet excluded by experiment. A recent 
analysis of experimental data from r-decay suggests a(M z ) « 0.129 |4lj . If we increase the 
number of flavors in our calculation we encounter large numerical uncertainties and do not 
obtain convergence for Nf > 5. This might be a signal that in the range 3 < Nf < 5 the 
above discussed transition of the system to the chirally symmetric phase takes place. 

All employed vertex constructions allow for nontrivial solutions of the quark equation 
corresponding to dynamical chiral symmetry breaking. However, similar to the quenched 



case and in accordance with the results of refs. ;:20, 211 the bare vertex construction and the 



CP-type vertex with d = 5 generate much too small quark masses compared with typical 
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FIG. 8: Displayed are the ghost and gluon dressing function, Z and G, the running coupling a, the 
quark mass function M and the inverse vector self energy 1/A. The calculations are done quenched 
and unquenched with Nf = 3 quarks in the chiral limit. The parameter d in the vertices is set to d = 0. 

phenomenological values. For d = we obtain good results for the quark mass, the pion 
decay constant and the chiral condensate, whereas the choice d = 0.1 leads to somewhat 
large values. It is interesting to note that d = of all values is preferred as in this case the 
quark equation resembles most the fermion equation of QED. 

In Fig. |Hlwe display the ghost, gluon and quark dressing functions corresponding to the 
unquenched and quenched cases with d = from table |H] |53j]. We find different anomalous 
dimensions in the ultraviolet corresponding to the change from Nf = to Nf = 3, c.f. 
eqs. (|50|) . ([SB)). ()89p . As expected from the infrared analysis of the quark-loop the back- 
reaction of dynamical quarks in the gluon equation does not affect the infrared behavior of 
the ghost and gluon dressing functions. Consequently the infrared fixed point of the running 
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x=p 2 [GeV 2 ] x=p 2 [GeV 2 ] 

FIG. 9: Results from the unquenched calculation with Nf = 3 massive quarks. We used the renor- 
malized masses m u / d (2GeV) 4.6 MeV and m s (2GeV) rj 96 MeV. Again we chose d = for both 
vertices. 

coupling is the same as in pure Yang-Mills theory. 

Our results for the case of explicitly broken chiral symmetry are shown in Fig. We 
choose Nf = 3 with renormalized quark masses corresponding to m M /^(2 GeV) = 4.6 MeV 
and m s (2GeV) = 96 MeV within our renormalization scheme. These masses are well in 
the range suggested by the Particle Data Group |34| , however, they should not be identified 
directly as the PDG employs an MS'-scheme. Compared to the chiral case the behavior 
of the ghost and gluon dressing functions hardly changes. For the quark mass function we 
obtain the irregular asymptotic solution in the ultraviolet as expected. 

For further use, e.g. in phenomenological calculations, we provide fits to our results for 
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1.22 


1.00 


1.33 


2.01 
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1.83 


1.23 














96 


98 


528 


2.96 


1.03 





A-\Q,s) 


hi 


h 2 


h 3 


h,4 




0.638 


0.515 


0.00688 


0.562 


0.00681 




0.671 


0.302 


0.00139 


0.318 


0.00137 



TABLE III: Parameters for the fits to the unquenched results with N f = 3, 5 = -0.25, j m = 12/27 
and fio = 27/3, using the CP-vertex with d = 0. The renormalization point s = 497GeV 2 is determined 
by the condition a(s) = 0.2. Note the change in Aqqq as compared to the case of chiral quarks. 



the quark propagator employing the fit functions 



M(x) 



9i + (x/Aqcd] 



92 



9iM(0) + 



m 



ln(x/A 2 



QCD) 



{x/A 2 



QCD) 



7m 



[x/A 



QCD) 



lA( r i [A(0, s)] + hi (x/A 2 QCD ) + h 2 (x/A 2 QCD Y 
[ [X,S)l ' l + h(x/A 2 QCD ) + h A (x/A 2 QCD ) 2 



(94) 



(95) 



with x = p 2 and the six parameters 91,92, h%, h 2 , h 3 , h^. We used the renormalization point 
independent current-quark mass rh, which is related to the renormalized mass M(s) by 



rh = M{s) 



Hs/^qcd) 



7m 



(96) 



to one loop order. For the running coupling, the ghost and the gluon dressing function we 
use the form 'Fit B', given in eq. ()38|) and the fit functions from eqs. (141 j) . In table ILTTl we 
give our values for all parameters as well as the numerical results for M(0) and [A(0)] _1 . 
Note that the scale Aq§^ is different to the corresponding scale in the chiral limit due to 
the different ultraviolet behavior of the quark-loop when quarks with non-vanishing bare 
masses are employed. When plotted the fits are virtually indistinguishable from our results 
in Fig. El 

Unquenched lattice calculations employing dynamical quarks are complex and time con- 
suming (43). To our knowledge such simulations for the propagators of QCD have not yet 
been performed. From our results in the Dyson- Schwinger approach we do not expect drastic 
differences between quenched and unquenched propagators on the lattice. Note, however, 
that our calculation includes quark-loop corrections to the gluon self energy but not higher 
order vertex corrections like mesonic loops. In a model calculation of the pion charge radius 
such loops have been estimated to contribute roughly at the order of ten percent 3 ■ 

Finally, we investigate possible positivity violations in the gluon and quark propagators. 
According to the Osterwalder-Schrader axiom of reflection positivity a two-point correlation 
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FIG. 10: Here we display the one dimensional Fourier transforms of the gluon propagator, D^t^p 1 ), 
and the scalar and vector parts of the quark propagator, crs{—t,jp) and ay(— t^p 2 ). For the three 
momentum we chose p 2 = 0. The time t is given in internal units. We observe violation of reflection 
positivity for the gluon propagator but not for the quark propagator. 



function S of Euclidean field theory has to fulfill the condition |4 

dtdt' f(t',p) S(-(t + t'),p) f(t,p) >0 



(97) 



if a physical particle is described. Here / are complex valued test functions. A violation of 
this condition signals the absence of the corresponding particle from the physical spectrum 
of the theory, i.e. the particle is confined. The one-dimensional Fourier transform S(t,p) of 
the propagator S(po,p) is given by 



S(t,p) := 



dpo 
2tt 



S(,Po,p)e 



ip t 



(9f 



Provided there is a region around t where S(—to,p) < one can choose a real test function 
f(t) which peaks strongly at to to show positivity violation. In the following we chose p — 0. 
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In the first diagram of Fig. we display the Fourier transform of the nontrivial part 
D(p 2 ) = Z(p 2 )/p 2 of the gluon propagator. Clearly one observes negative values on a large 
interval. The resulting positivity violation for the transverse gluon propagator in Landau 
gauge is a clear signal for gluon confinement. This corroborates previous findings in the 
quenched approximation 0. These positivity violations have also been observed in lattice 



studies, see ref. jl2| for recent corresponding results or ref. |46( for a review. 



In the quark propagator positivity violations have been found in model studies, which 
solve the quark Dyson- Schwinger equation with an ansatz for the gluon propagator as in- 



put (see |2. 47, 48] and references therein). Similar violations have been found in (2+1)- 
dimensional QED ji^. For the model of ref. |33[ or a propagator with complex conju- 
gate poles we tested our numerical routines and found positivity violations very easily. 
Contrary to these findings we do not observe positivity violations for the quark propaga- 
tor from the coupled set of DSEs. The lower panel of Fig. El shows our results for the 
Fourier transform of the vector part <Jy{p 2 ) = A(p 2 ) / (p 2 A(p 2 ) + B(p 2 )) and the scalar part 
^sip 2 ) — A(p 2 ) / (p 2 A(p 2 ) + B(p 2 )) of our solutions for the quark propagator employing four 
different quark-gluon vertices. All our solutions appear to be positive definite at the present 
level of numerical accuray. However, a more accurate study is required to settle this point 
jo"ll ]. Furthermore note that even if confirmed our findings are not in contradiction with the 
absence of quarks from the physical spectrum of QCD as violation of positivity is a sufficient 
but not a necessary condition for confinement. 



SUMMARY 

We have presented solutions of the (truncated) Dyson-Schwinger equations for the prop- 
agators of Landau gauge QCD. We first concentrated on the Dyson-Schwinger equation for 
the quark propagator. We proposed several ansdtze for the quark-gluon vertex which consist 
of an Abelian part carrying the tensor structure of the vertex and a non-Abelian multiplica- 
tive correction. Our guiding principles for the construction of these vertices have been two 
important conditions on the truncated quark equation: it should be multiplicatively renor- 
malizable and recover perturbation theory for large external momenta. In our truncation 
scheme the quark mass function is, as required from general arguments, independent of the 
renormalization point and has the correct asymptotic behavior for large momenta. 

In the quark equation both the ghost and gluon dressing function show up at least implic- 
itly. In quenched approximation, which is suitable to compare to lattice results, we employ 
solutions of the ghost and gluon Dyson-Schwinger equations taken from ref. Q. In a second 
step we included the back-reaction of the quarks on the ghost and gluon system and solved 
the quark, gluon and ghost Dyson-Schwinger equations self-consistently. 

All our solutions exhibit dynamical chiral symmetry breaking. However, only carefully 
constructed vertex ansdtze have been able to generate masses in the typical phenomeno- 
logical range of 300 — 400 MeV. Constructions with an Abelian part satisfying the Abelian 
Ward-identity are superior to other vertex ansdtze. We obtained very good results for the 
quark mass, the pion decay constant and the chiral condensate by employing a generalized 
Curtis-Pennington 123] vertex. In the chiral limit both, the quark mass function and the 
vector self energy coincide with recently obtained lattice results [l3l . I3 | within the numerical 



uncertainty. This agreement confirms the quality of our truncation and in turn shows that 
chiral extrapolation on the lattice works well. 

In the unquenched case including the quark-loop in the gluon equation with Nf = 3 
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light quarks we obtain only small corrections compared to the quenched calculations. In 
particular for the case of dynamically generated quark masses the quark-loop turns out to 
be suppressed in the gluon equation for small momenta. We thus showed on the level of 
our truncation that the Kugo-Ojima confinement criterion 3^, and Zwanziger's horizon 
condition [1, H(| are satisfied in Landau gauge QCD. 

Furthermore we searched for positivity violations in the gluon and quark propagators. 
We confirmed previous findings |5| that the gluon propagator shows violation of reflection 
positivity. Thus the gluon is not contained in the physical state space of QCD. We did not 
find similar violations for the quark propagator. This issue is currently investigated in a 
more detailed study p5ll ]. 
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APPENDIX A: ANGULAR INTEGRALS 

The angular integrals employed in the ultraviolet analysis are given by 
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where the squared momentum z is defined as z = (p — q) 2 = x + y — 2Jx~y cos(9) 
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For a discussion of the relation of a MOM-scheme to the MS'-scheme see section IV of ref. |2f 
Based on a three loop calculation the authors obtained the relation A-p^ ~ 0.346A A ^ M . In 
our case this would result in the value A-p^ = 0.246 which for zero flavors is even on the low 
side. 

[53] Graphs for the choices d = 5 and d = 0.1 can be found in ref. [421 ] . 
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